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ABSTRACT 

The orbits about Lagrangian equilibrium points are important for scientific investigations. 
Since, a number of space missions have been completed and some are being proposed by various 
space agencies. In light of this, we consider a more realistic model in which a disk, with power-law 
density profile, is rotating around the common center of mass of the system. Then, we analyze 
the periodic motion in the neighborhood of Lagrangian equilibrium points for the value of mass 
parameter < /i < i. Periodic orbits of the infinitesimal mass in the vicinity of equilibrium are 
studied analytically and numerically. In spite of the periodic orbits, we have found some other 
kind of orbits like hyperbolic, asymptotic etc. The effect of radiation factor as well as oblateness 
coefficients on the motion of infinitesimal mass in the neighborhood of equilibrium points are 
also examined. The stability criteria of the orbits examined by the help of Poincare surfaces of 
section (PSS) and found that stability regions depend on the Jacobi constant as well as other 
parameters. 

Subject headings: Periodic orbits: Photogravitational: Oblateness: Disk: Chermnykh-like problem: 
Poincare surfaces of section. 



1. Introduction 

For the last few years, many authors have stud- 
ied the periodic orbits of restricted three body 
problem (RTBP) due to its wide range of ap- 
plications in sp ace dynamics. One of them was 
Plummeil (|l90l f) who studied the planar circular 
RTBP with arbitrary mass parameter /i and found 
that there are two families of pe riodic motion near 
the Lagrangian points. It was ISzebehelvl (jl96lT ) 
who described complete results regarding the peri- 
odic motion of planar circular RTBP. A systematic 
classification of periodic orbits for < /i < ^ in 
th e neighbo r hood of Lagrangian points presented 
bv lBrouckd (Il968h . The periodic orbit about tri- 
angular point with ma ss parameter as critic a l mass 
of Rout h studied b y iMever and Schmidtl (|l97ll ) 
whereas iMarkellos ( 1974) iiiyestiga ted the prob- 
lem numerically. iRagos et all (|l99l[ ) examined the 
periodic motion around the collinear equilibrium 
points in the photogravitational RTBP. However, 



lElipe and Laral ([19971) discussed the same motion 
by takin g both primar ies of RTBP as radiating. 
Further, iPerdiosI (|2003l ) studied critical symmet- 
ric periodic orbits of RTBP t aking one primary as 
an oblate spheroid. Whereas, iHenon (l2003l l2005l) 
investigated new families of periodic orbits during 
the study of Hil l's problem of the t hree body. On 
the other hand, iMittal et al (l2009l) examined the 
periodic motion of RTBP with oblateness by tak- 
ing displacements along tangent and normal to the 
mobile coordinates. 

To analyze stability of the orbits, there are 
several methods like Lyapunov's , Poincare Ma p 
or Poincare surfaces of section (jPoincard Il892r) , 
Henon's horizontal-vertical indices etc. The 
PSS methods has been used by many authors 
like (iRagos et a]||l997t IWinter and Murravl [l997 : 



Winter 200(D , to study the nature of orbits around 



the equili brium points in the restric t ed th ree body 
problem. ISafiya Beevi and Sharma ( 2011 ) studied 
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the periodic orbits in the Saturn-Titan problem 
using the numerical technique of PSS and found 
that the orbits around Saturn remain around it 
and their stability increases with the increase in 
the value of Jacobi constant. 

In the present paper, we have analyzed the peri- 
odic motion in the generalized photogravitational 
Chermnykh-like problem with power law density 
profile of disk which is rotating around the com- 
mon center of mass of the system with mass ratio 
< /i < ^- The Chermnykh-like problem has a 
number of applications in different areas such as 
celestial mechanics, chemistry, extra solar plan- 
etary system etc ( Gozdziewski and Maciejewskjl 



1999 ; iStrand and Reinhar dt 1979 ; Rivera and Lissauei 
2000: I Jiang and id l200ll). This p r oblem was 
first time studied by Chermnvkh' ("1987 1) lat- 
ter ma n y authors (Gozdziewski .1998.: .Paoadakis 



2005allbl: Jiang and Yehll200fiHYeh and Jianel200fil: 



Kushvahll2008l) have anaysed the different aspects 
of the problem theoretically as well as numeri- 
cal ly. The propo s ed pr oblem has been studied 
by kushvah eta] (|2012l ) in the consequence of 
existence of equilibrium points and their linear 
stability. 

The numerical computati on has been don e with 
the help of Mathematica ® IWolframI (l2003l) soft- 
ware at necessary places by taking initial dis- 
placements ^0 = riQ — 0.001 and initial velocities 
Co = Vo = 0.01. The parametric values used in 
all numerical results are qi = 0.75; A2 — 0.0025; 
c = 1910.86; a = 1; 6 = 1.5 and h = lO"* unless 
otherwise stated. To find the PSS of the prob- 
lem, we have used Event Locator method which 
is built into Mathematica® and works effectively 
as a controller. Since, the errors in a numerical 
method, generally, depend on step-size used which 
is controlled internally. Therefore, in a lower or- 
der integration routine error will be greater and 
hence more inaccuracy in results. As, a numerical 
method yields approximate results, consequently, 
there is an unavoidable error called as round off 
error which is the difference between approximate 
result and exact result. This error occurs due to fi- 
nite precision representation of floating point num- 
bers in the computer arithmetic. Since, we do not 
know the exact solution of the problem so, we can 
only control the errors by increasing working preci- 
sion. On the other hand, a high working precision 
is unfavorable to computer and takes more time 



for execution of the results. Therefore, for a large 
time t round off error will be large. 

The work is presented in six sections in which 
section(I2]) contains the mathematical formulation 
of the problem and equations of motion. Sec- 
tion ([3]) includes the variational equations in the 
vicinity of equilibrium points. Section(|4]) and ([5]) 
describe the periodic motions of the massless body 
in the neighborhood of collinear and triangular 
equilibrium points respectively. In section(|ni), we 
have discussed the stability criteria with the help 
of PSS technique while section([7]) concludes the 
paper. 

2. Mathematical Formulation of the Prob- 
lem and Equations of Motion 

It is supposed that mi and m2{mi > TO2) be 
the masses of first (radiating body) and second ( 
oblate body) primaries respectively; Md be the to- 
tal mass of the disk, rotating around the common 
center of mass of the system. The power law den- 
sity profile of the disk having thickness h w 10"'' 
is p(r) — where p is natural number(here we 
take p = 3) and c is a constant determined with 
the help of mass Md- The forces, governing the 
motion of infinitesimal body m are gravitational 
attractions due to both primaries as well as due to 
the disk. The radiation pressure force and oblate- 
ness coefficient are also taking into account. It 
is also assumed that the effect of the infinitesi- 
mal body on the motion of the remaining system 
is negligible. The unit of mass is taken in such a 
way that G{mi-\-m2) = 1; unit of distance is taken 
as the separation between the primaries whereas 
unit of time be the time period of the rotating 
frame Oxyz. The common center of mass of the 
primaries is taken as the origin O of the rotating 
frame which is fixed relative to the inertial sys- 
tem. Let, P{x, y, 0), A{-n, 0, 0) and B{1 ~ 0, 0) 
be the co-ordinates of infinitesimal body, first pri- 
mary and second primary respectively, with re- 
spect to the rotating frame, where /i = m^m2 
mass parameter. 

Therefore, the equations of motion of the 
i nfinitesimal body in xy-plane is written as 
(|Kushvah et alll2012l) : 
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The potential function 51 is given as: 

r2 



n = 



2 ri 
where, y is potential of the disk; r 
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(1) 
(2) 



(3) 



mass reduction factor q\ = (1 — -p-), where Fp 
and Fg are the radiation pressure and gravita- 
tional attraction force of the radiating body re- 

specti yely; A2 = \m " is the oblateness coeffi- 
cient ( McCuskevlll963l ). where Rg and Rp are the 
equatorial and polar radii of the same body re- 
spectively and R is the distance between the pri- 
maries. The mean motion n of the system is given 

by ^Jqi + 1^2 - 2/f,(r), where /^(r) is the gravi- 
tational force due to the disk. 

The potential as well as gra vitational force of 
the disk are given as (jjiang and Yeh.2006h : 



y = -4 



F{Op{r'y 



dr' 
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p{r')r' 



(4) 



dr', (5) 



where F{() and E{Q are elliptic integrals of first 
and second kind respectively, r' is the disk's refer- 
ence radius and C — 2^^. Now, with the help of 
expansion of the elliptic integrals for a < r' < b in 
dU and ([5]) and then selecting appropriate terms 
relative to r, the simplified form of fb{r) is given 
as: 



Mr) 



-TTc/l 



2{b-a) , 3(log^) 



abr^ 



8r3 



(6) 



where, a and b are inner and outer radii of the disk 
respectively. It is assumed that the gravitational 
force fb{r) is radially symmetric and hence ffb{r) 
and ^fb{r) are taken as the components of the 
force fbir) along x and y axes respectively. From 
equation ([!]) and ^ , we obtain the Jacobi integral 
of the problem given as: 



C = -x^ - f 



2U, 



(7) 



3. Variational Equations in the Neighbor- 
hood of Equilibrium Points 

To analyze the possible motions of the re- 
stricted body in a small neighborhood of the 
equilibrium points {xe,ye), we first make an in- 
finitesimal change ^ and 77 in its coordinates i.e. 
X = Xe + ^,y = ye + V such that the displacements 



^ = Pe^\ r, = Qe 



At 



(8) 



are very small, where P, Q and A are parame- 
ters to be determined. Substituting these coor- 
dinates into equations dU) and ([2]), we get two 
differential equations of second order in vari- 
able £ and 77 known as variational equations 
(jMurrav and Dermottll2000[ ): 



2nf] : 
-2< 



(9) 



where superfix indicates the values are com- 
puted at the equilibrium point. Again, putting 
£ = Pe^*, 77 = Qe'*'* into the equation Q and sim- 
plifying them, we have 



{X'-nlJP+i-2nX^n°^)Q^0, 



(10) 

(2nA-f7°jP+(A2-flO^)Q = 0. (11) 

Now, the condition of nontrivial solution is that 
the determinant of the coefficients matrix of above 
system should be zero i.e. 



A2- 
2nA 



XX 

)0 



-2nA 



= 0. 



So, from above we obtain a quadratic equation in 
A^ known as characteristic equation: 



A^ + (4n2 - - nl)X' + 



XX 

2 



i^xx^VV ^ xy) ~ 0- 



(12) 



The four roots of characteristic equation ([T^ play 
a crucial role to determine the form of orbits in the 
basin of equilibrium points. Now, the second order 
derivatives of potential function Q with respect to 



where C is Jacobi constant. 
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X and y are given as: 



where, 



^7J 



3(a: + At)2 



3^2 

2ri 



TTch 



2(6 -a) 
ah 



r 

TTchx'^ 



3/i(a; + — 1)'^ 



6(6 -g) 7(log^) 
a6 r 



5^2 

2ri 



(13) 
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3^2 

2ri 



3y^ 



TTch 



2(6 -a) 



a6 



7(log^) 



3y^ 

To 



5^2 

2rl 



-Kchy^ 



6(6 -a) 7(log^) 



a6 



,(14) 



3gi(l - ^i){x + fj.)y 



3n{x + n - l)y 



5^2 

2rj 



TTchxy 



6(6 -g) 7(log^) 
a6 r 



(15) 



gi(l - /^) M 

3fj,A2 2'Kch{b — a) 



2\{XL^+IX-1)\^ 

77rcft(log^) 
4(:el.)4 • 



ab\xLi 



Now, the last term of the characteristic equation 
are computed numerically by using the para- 
metric values stated earlier and results are plotted 
in figure (H]). It is found that 



(17) 



for all the collinear equilibrium points XLi,i = 
1,2,3 under the condition of mass parameter /i 
w hereas the abov e inequality holds for all values of 
jjL ( Brouckell 19681 ). In other words in case of Li and 



L2 above inequality holds good for 0.01 < /i < i 
while in case of L3 it is true for 0.13 < < i. 
Since, the zero degree term in the characteristic 



\ 



\ 



Fig. 1.- Values of = ^^x^l^ ~ [n^f 
with respect to mass parameter /i: {ajx^-^ — 
0.813609; (6)xl, = 1.05667 and {c)xl^ = 
-0.823420 



4. Periodic Orbits in the Neighborhood of 
the CoUinear Points 

The values of second order derivatives of VI at 
the collinear equilibrium points (x^. , 0), i — 1,2,3 
are given as: 



^''xy 



3B 



{xL, + M - 1)1' 



0, 



(16) 



equation ([T^ has always a negative sign for these 
values of /x. Hence, inequality pT|) guarantees that 
we always have two real roots of opposite sign in 
under the restricted value of /i. That is 



a 



(18) 



Thus, the four roots Aj,j — 1,2,3,4 of the char- 
acteristic equation are given as 



Ai 
A3 



-a, A2 = 
A4 



-z/3. 



(19) 



Hence, from equation (|5]) it is clear that the or- 
bits corresponding to the roots Ai and A2 are of 
the exponential type. But orbits corresponding to 
roots A3 and A4 are periodic with period 
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Now, the most general solution (which is of the 
complex type) of the variational equations ([9]) cor- 
responding to the roots ([T9l) can be written as: 



^ - Pie"* 



P2e-"' + Pae'-^* + P 



46 



(20) 



where Pj,j = 1,2,3,4 are four arbitrary integra- 
tion constants, whereas constants Qj,j — 1, 2, 3, 4 
are related to previous constants by the linear 
equation (fTO|) . That is 



2\i 



P, J = 1,2,3,4. 



(21) 



Let us suppose that 



(22) 



where 

7 = 



1 

2^ 



7^c/i(logf) 
4(xlJ4 



- 3B 

fc^ 



3m^2 



|(xi. +M-1)P 



(23) 



1 

2^ 



3P 



3AtA2 



(xi. +M-1)|^ 



7^c/i(logf) 



(24) 



4(xlJ4 

Then, the general solution ([20)1 takes the form 

C - Pie"* + P2e-"*-HP3e*^* + P4e~*'3*, 
= 7(Ae"*-P2e-"*) 

+i6 (Pae'^* - P^e-'f^*) , (25) 

Also, the most general real solution of variational 
equations ([9]) can be written as: 

£, = Pi cosh at + P2 sinh at + P3 cos (3t + 
P4 sin /3t, 

?7 = 7 (Pi sinh at + P2 cosh ai) — 6 (P3 sin /3i 
-hP4 cos I3t) , (26) 

where Pj,j = 1,2,3,4 are arbitrary real integra- 
tion constants. Now, taking appropriate initial 
conditions for the perturbed motion in the basin 
of coUinear point, we have different type of orbits. 



4.1. Periodic Orbits 

It is supposed that initial displacements and ve- 
locities are Vo^ Co and rjQ respectively. Then at 
t = 0, we have from the equations ([20]) : 

eo = Pl+P2 + P3+P4, 
??0-7(^l-^2)+i<5(P3-P4), 

eo = a(Pi-P2)+*/3(P3-P4), 

rio - aj (Pi + P2) + 1(35 (P3 + P4) . (27) 

Again, from the equation (f25|) it is clear that for 
periodic motion. Pi and P2 must be zero i.e. Pi = 
P2 = . Thus, from first two equations of (|27p . 
we get 



P3 



26 



1 ^4 — 



26 



(28) 



Using above, the equation p5|) reduces to follow- 
ing form: 



6^,0 - irjo 



i6 



25 

6^0 - im 



25 



-ipt 



25 
26 



JI3t 



-iPt 



(29) 



Now, simplifying these equation with the help of 
Euler's relations, we obtain 



C = Co cos /3t + ^ sin pt, 
6 

?7 = 770 cos pt — (5^0 sin [3t. 



(30) 



The parametric equations (j30p represent the peri- 
odic orbits of the infinitesimal mass in the neigh- 
borhood of coUinear points Li with periods Ti = 
i = 1,2,3 which are depicted in figure ([3]). 
Now, eliminating parameter t from these two equa- 
tions, we get 



= 1. 



(31) 



which is the equation of an ellipse with center at 
coUinear points Li, L2 or L3 and axes parallel to x 
and y axis of the rotating frame. Also, the ellipse is 



bounded by rectangular region ^ = iy (^q -I- |8-) 
and rj = ±^/(WQ+r|^, where 6 defines the shape 



5 



of the ellipse. The eccentricity of this ellipse (jST 
is given as: 



(32) 



In figure frame (a) shows that the eccentric- 





Fig. 2. — Variation of eccentricity 

ity of ellipse lie in between interval (0, 1) and vary 
with the value of 6. In other words, ellipticity 
of the orbit get increases by increasing the value 
of 5 with condition 6^ > 1 whereas it decreases 
when 6^ < 1. Accordingly, the axes of the or- 
bits get change with time. From equation it 
clear that S is the function of parameters qi , A2 
and disk's radii and hence, the shape of orbit de- 
pends significantly on these parameters. Frame 
(6), shows that variation of eccentricity of the hy- 
perbolic orbit (p8| with the value of 5. 

Alternatively, if we take Pi = P2 = in the 
real solution (pS)) . we have 

^ = P3 cos I3t -f P4 sin pt, 

f^^-S (P3 sin /3t + P4 cos /3t) , (33) 

which also represents a periodic motion. Simplify- 
ing equation ([55]) by taking displacements £, = S,o, 
77 = 770 at i = 0, we get 



C = $0 cos/3i; 7] = sin/3i; 



(34) 



which is a one parameter family of an ellipse hav- 
ing center at coUinear points Li, L2 or L3 and in- 
tersecting ^-axis perpendicularly. Due to the pos- 
itive value of constant S motion on these ellipse 
will be retrograde. 

In figure the frames (a), (5) and (c) in- 
dicate that periodic orbits in the neighborhood 
of coUinear points Li, L2 and L3 with periods 
1.25876, 0.843075 and 3.76437 which are drawn 
at mass parameter fj, = 0.05, 0.05 and 0.15 re- 
spectively. It is seen that orbits get a deflection 



from their original path when we change the val- 
ues of parameter qi or ^2- In other words, the 
effects of radiation pressure and oblateness on the 
motion of infinitesimal mass are considerable. It is 
found that, generally the nature of the motion is 
unaffected but there is a miner change in time pe- 
riod of the orbits shown in Table [T] In this table, 
last three columns indicate the time periods of the 
periodic orbits in the vicinity of coUinear points 
Li, L2 and L3 respectively. It is clear that the 
time periods of the orbit decrease with increasing 
the value of qi or A2. In other words, time periods 
of the orbits in the basin of coUinear points de- 
crease with an increase in radiation pressure. The 
oblateness coefficient also reduces the time periods 
of orbits. That is, if A2 increases then the motion 
around equilibrium points becomes fast. 




Fig. 3. — Periodic orbits in the neighbor- 
hood of coUinear equilibrium points: {a)xL-^^ = 
0.813609; (6)xl, = 1.05667 and {c)xl., = 
-0.823420 



Table 1: Effect of radiation factor and oblateness 
on the motion of infinitesimal mass 
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A2 


= 0.0025 




0.75 


1.25876 


0.843075 


3.76437 


0.85 


1.25553 


0.842186 


3.53663 


0.95 


1.25230 


0.841297 


3.34299 


A2 




qi = 0.75 




0.0020 


1.25878 


0.843082 


3.76995 


0.0025 


1.25876 


0.843075 


3.76437 


0.0035 


1.25871 


0.843061 


3.76065 


0.0040 


1.25869 


0.843055 


3.75881 
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4.2. Hyperbolic Orbits 

Now, to discuss the unbounded motion of in- 
finitesimal mass, we assume that the arbitrary 
constants P3 and P4 in equation (j25p are zero. 
Again, taking the same initial displacements and 
velocities as in case of periodic orbits. Then, from 
first two equations of ([?f)) . we find the following 



Pi 



7C0 + Vo 
27 



7C0 - Vo 
27 



(35) 



Substituting these values of arbitrary constants in 
the equation (PO]) . we get 



76 + "iyo 
27 

7^0 + m 
27 

76 - im 

2a 



7C0 - iyo 

27 



(36) 



Euler's relations reduces the equation in fol- 
lowing form 

^ = ^0 cosh at + — sinh at, 

7 

1] = 7C0 sinh at + rjo cosh at. (37) 



The parametric equations p7p represent path of 
the infinitesimal body in the vicinity of collinear 
points for the appropriate values of fj, as described 
earlier. Eliminating parameter t, we have 



V 



= 1, 



(38) 



which is the equation of hyperbola with center 
at Li, L2 or L3 and axes parallel to the coor- 
dinates axes of the rotating frame. The hyper- 
bola ((38)) is bounded by the rectangular region 



e = ±^ieo - ^) and 77 = ±x/(7^eo'+'7o), where 
7 defines the shape of the hyperbola. The eccen- 
tricity of this hyperbola is given as: 



(1 
(1 



^)^; 7'>i, 



-7 



< 72 < 1. 



(39) 



and P4 in the equation ([261) as equal to zero. If 
we changing the time origin then motion can be 
represented as: 



^ = ^0 cosh at; -q — 7^0 sinh at; 



(40) 



which represent the one parameter family of hy- 
perbolas with center at Li, L2 or L3 and crossing 
the <^-axis at right angle. 

4.3. Asymptotic Orbits 

The asymptotic motion is possible for the par- 
ticular values of the real arbitrary constants in 
the equation (PB)) . In other words, if we choose 
P3 = P4 = and Pi = ±P2 then the solution is 
written as: 



C = ±Coe 



±at. 



77 = ±7^0 e 



±at. 



(41) 



which represent a rectilinear path of the infinitesi- 
mal body in the basin of the collinear points. This 
path may be considered as a hyperbola degener- 
ated in its asymptotes. Plus minus sign in the 
equation (|4T|) , indicate that the infinitesimal body 
going away from or coming towards one of the 
collinear points Li,i — 1,2,3 asymptotically with 
time t. For each collinear equilibrium point, there 
are four possible asymptotic orbits according to 
the signs of equation (|4ip . Again, due to sym- 
metry of the RTBP relative to x-axis, only two 
asymptotic orbits are really fundamental while 
other two are images of them. In other words 



(42) 
(43) 



are two fundamental out going orbits, while two 
are incoming orbits 

? = eoe-"*, ,7 = -7^06-"*, (44) 
C = -Coe-"*, V^lCoe-"' (45) 

are images of (|^ and (^5]) respectively. 

5. Periodic Orbits in the Neighborhood of 
Triangular Equilibrium Points 

The coor dinates of the triang ular points £4^5 
are given as (jKushvah et a]||2012i ): 



Alternatively, the hyperbolic motion in the 
vicinity of a collinear point can be examined by 
choosing the values of real arbitrary constants P3 



X = ^ /i 



,y = ±ql [1 - 2i_ + (2 - ql)ai + a^]^ 



[qlai - (72) 

2 

i? 
4 



(46) 
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where ui and a2 are given as: 



2TTch{b — a) 



3^c/i(logf) 



8{/.2+qf(l-^)}2 
1 



a'2 



3(1 + 1^2 



1 — n H 



2TTch(b — a) 



- + 

ab{^l^ + qlil^^l)}-2 

37refe(logf) 

8{M2 + gf(l-M)PJ 

Now, using parametric values mentioned earlier in 
addition with triangular point P5|) into the equa- 
tions ([T^ , (ITil) and ([T5|) , we have found the values 
of r2°^, riyj^ and V,^y. Again, substituting these 
values in the equation we get another equa- 
tion in the variable A and parameter /i which has 
a very complicated form. Again, expanding this 
equation by Taylor's series about the mass param- 
eter /i = and taking up to the first order term, we 
found the following new characteristic equation: 



+ (3.24710 + lM976fi)X^ 
7.38874^4- 2.36659 = 0. 



(47) 



Since, the nature of roots of the quadratic equation 
(|47)) in A^ depends on the sign of its discriminant 
D which is given as: 



D = (3.24710-1- 1.84976^)- 
4 (7.38874// -f 2.36659) . 



(48) 



Thus, we have the following three cases. 

5.1. When D = 0; /i = /i^ 

This value of the discriminant directly corre- 
sponds to the value of critical mass fi^ of this prob- 
lem. In other words 

3.42160 {^i - 5.06476) (/i - 0.062165) = 0. (49) 

So, we have the critical mass fic — 0.062165 be- 
cause fi = 5.06476 is not possible as < fi < \. 
For this value of /i, the characteristic equation (|T7)) 
has four imaginary roots Ai_2,3,4 of equal modulus 
given as: 



ai: 



-ai. 



(50) 



where a = 1.29655. Hence, the solutions of equa- 
tion (lU corresponding to distinct roots will be 
bounded however due to repeated roots, the gen- 
eral solution of the variational equation ^ will 
be unbounded. The most general solution of vari- 
ational equation ^ can be written as: 



^={Pi+ P2t) e'"* 
= (gi + Q2t) e'"* - 



{Q3 + Q4t)e''"\ 



(51) 



which is of complex type, where Pj,j = 1,2,3,4 
are arbitrary constants whereas Qj,j = 1,2,3,4 
are related to them by the equation (fT0|) . That is 

/a? - nt, 



2A, - n% 



P„j = 1,2,3,4. (52) 



The most general real solution of variational equa- 
tion ([9]) can be written as: 

^ = P cos {at + (p) + P't sin [at + cj)'), 

T] ^ Q cos {at + tp) + Q't sin {at + ^p'), (53) 

where P' , Q' , (j)' , ip' are related to the four arbi- 
trary integration constants P,Q,(j),Tp respectively, 
which are given as: 



P 



Pf 



Pi Q 



Ql + Q 



~ tan 



^ ) , -0 = tan 



-I'gi 1.(54) 



Now, to discuss the periodic motion, we have taken 
P2 = P4 = 0. Hence, from the equation ([5T|) we 
have 



^ = Pie*"* -l-Pae"'"*, 
7/ = 7iAe"* - 73^36-'"*, 



(55) 



where 7i = p|- and 73 = Again, if ^0 and 770 
are the initial displacements then equation (j51|) at 
t = provides: 



Pi 



736 - '^0 . 
73-71 



P3 = . (56) 

71 - 73 



The equation (j55p represents the periodic orbits 



with period , „f?L^.. of the infinitesimal mass in the 

^ 1.29655 

vicinity of triangular equilibrium £4 which are de- 
picted in figure (jlj. 

Figures (U) and ([S]) shows that the orbits of the 
infinitesimal mass plotted at different time. From 
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these figures, it can be seen that the orbit have a 
regular elliptic shape and retaining its path for a 
long time t > i.e. ii t = SOOT, then the shape of 
the orbits remains same as for t = 5T. It is also 
seen that if we draw the parametric equation (|5ip 
by assuming Pi — P3 ^ 0, then the path of the 
infinitesimal mass is moving spirally outward with 
time. 

5.2. When D > 0; /ie(0,/ic) 

In this case, due to positive sign of D, the char- 
acteristic equation (|47p has four imaginary roots 

Ai^2 = A3, 4 = ±/3i, (57) 

where a and /? are positive real numbers whose 
values depend upon the mass parameter fi. Due 
to two different imaginary roots ai and (3i, there 
are two periodic motions of the restricted body 
with periods ^ and The general complex so- 
lution of the variational equation ^ in this case 
is written as: 

^ = Pie'"* + Fae"*"* + Pae'^* + Pie^"^\ 
= gie-* + Qae-^"* + Q^e'^' + Q4e~'^\ (58) 

where Pj,j = 1,2,3,4 are four arbitrary integra- 
tion constants and Qj,j — 1,2,3,4 are related to 
previous constants by means of relation ([52|) . The 
more general real solution of ([9]) may be written 
as follows: 



^ = P cos {at + (f))+Q sin {l3t + ^p), 
V = |7i|^'cos (at + + ei) 

+ |73|Qsin(af + V-e3), (59) 




Fig. 4. — Periodic orbits in the neighborhood of 
U: (a) < < < T(= ^^) (b) < i < 5T 



where P,Q,(I),'>P ^-re four real arbitrary constants 
given as: 

P = \IPi + Ph Q = \IPi + Pl 

<^ = tan-i (^^), V' = tan-i (§)-(60) 

The constants |7i|, I73I and ei, £3 are modu- 
lus and amplitudes of the complex number 7^ = 
= 1,3 respectively. The equation (1551) rep- 
resents composed form of two periodic motions 
known as long periodic motion and short periodic 
motion with periods Tq, = 2£ and Tp = ^ re- 
spectively. The amplitudes of these type of mo- 
tions are obtained from the integration constants 
Pj,j = 1, 2, 3, 4. The looping nature of trajectory 
of the infinitesimal mass shown in figures is ob- 
tained from the two different type of motions that 
contribute to the perturbed orbit about the equi- 
librium point. From figure ©(a), the time period 
of the resulting motion is T = 30.4 (as in our case 
a and /3 obtained for < n — 0.00001 are not inte- 
gers). Again, in case of ^ = f i where s and k are 
integers, the time periods of the resulting orbits 
is given by the LCM of Ta and Tp. For example, 
time period of the orbits shown in figure ^ are 
2tt whereas in frame (a), ^ = ^ is equal to ^ and 
in frame (6) it is equal to i. Also, the shape of 
this resulting orbit is bounded and has two loops 
in former case (figure [7]) while in later case (figure 
[9Ua&6)) have one and two loops respectively. The 
loops in a quasi-periodic orbit, are very helpful to 
inspect the order of resonance. 

Now, to study the separate motion of the in- 
finitesimal body, we have considered the motion 



17(0 r,Ct) 




Fig. 5. — Periodic orbits in the neighborhood of 
L4: (a) < t < 200T (b) < i < 500T 
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((59ll with Q = 0. That is 



£, = P cos {at + 4>), 
T] — |7i|Pcos {at + 



(1 



(61) 



which represents a one parameter family of el- 
hpse. Now, ehminating parameter t by trigono- 
metric simpHfication in between both the equation 
of (16T1), we get 



171 1 



2^77 9 2 

- — rcosei — P sm ei. 
l7i| 



(62) 



The equation ([62]) represents an elhpse which is 
bounded by the region ^ = ±P and r] ~ ±|7i|P. 

In Figure ® which is plotted at /i^ = 0.00001, 
frames (a) and (6) show that the orbits of infinites- 
imal mass in the vicinity of L4 with periods ^ 05102 
and ^ 1q^'j2 correspond to roots \i,2 = ±1.05102i 
and A34 = ±1.46372i respectively. From this fig- 
ure, it is also clear that orbits are elliptic and 
bounded for a long time t > 0(as in figure [5]) . 
On the other hand, we have depicted the shape of 
composed form of orbits of two periodic motions in 
figures ([7]) and (|8]) at different time t when ^ / f . 
Figure © for the case when ^ = f , where s and 



k are integers i.e. in frame (a) 



= i and 



frame (&) it is |. From figure ©(a), it is clearly 
seen that there is a second order resonance in the 
quasi-periodic orbit ([58)) . 



5.3. When Pi < 0; /^e(^c,^) 

In this case the four roots of the characteristic 
equation (|T7|) have the form 

Xi = a + iP; X2 = a — ij3\ 

A3 = —a + ijS] A4 = —a — ij3, (63) 



/ 0.02 

/ 0.01 


L4 / 


/■0.04 -0.02 

-0.01 

\ -0.02 
-0.03 


0.02 o.o4 



f(t) 



(b) 


'ao3 




„/ 


0.02 






0.01 










-0.04 - 


0.02 


0.02 0.04 




-0.01 






-0.02 






-0.03 





Fig. 6. — Periodic orbits in the neighborhood of 
U: (a) < i < P„(= (b) < t < T^(= 

1.46372 




Fig. 7. — Orbits in the neighborhood of P4 when 
^ = j: and s, k arc not integers: (a) < t < T{= 
30.4) = 7.IT/3 (b) < < < 2r 




Fig. 8. — Orbits in the neighborhood of P4 when 
^ = -I and s, k are not integers: (a) < t < ST 
(b) < i < 30r 




Fig. 9. — Orbits in the neighborhood of P4 when 
^ = I and s, k are integers: (a) < t < 27r (b) 
< t < 27r 
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where a and /3 are positive real numbers depend- 
ing on the mass parameter fi. These four roots are 
responsible for the motions of infinitesimal body 
in the basin of triangular equilibrium points. The 
motion is periodic with period ^ however the am- 
plitude is an exponential function of parameter t 
and a. In this case the most general complex so- 
lution of the equation (|9]) may be written as: 



77 = Qie^i* -f Qae 



^36' 



-,A3t 



-,A4t 



(64) 



where constants Qj,j = 1,2,3,4 related to four 
arbitrary integration constants Pj, j = 1, 2, 3, 4 by 
mean of ((52l) . Also, the most general real solution 
of ([9]) have the form: 

f = e"* (Pi cos + P2 sin I3t) + 
e""* (P3COS/3i-|-P4sin^t), 

7? = e"* (Qi cos pt + Q2 sin (3t) + 

e""*(Q3COS^t-h04sin/3i), (65) 

where real constants Qj, j = 1, 2, 3, 4 are the func- 
tion of arbitrary real constants Pj,j = 1,2,3,4. 
Now, if we take P3 = P4 = in equation (|65)) then 
we found that orbit is of asymptotic type which 
moves away from the triangular equilibrium point 
with time t. The period and amplitude of this 
orbit are ^ and e"* respectively. On the other 
hand, if we assume P^ = P2 = in equation (|55|) 
then the motion of infinitesimal body is similar 
to previous one provided amplitude e~"* should 
decrease with the time. 

In figure (fTU)) which is plotted at /ic = 0.0651, 
frame (a) and (b) show that the path of infinites- 
imal mass is of asymptotic type. It is also seen 
that the orbit is going away from center L4 with 



time t. Figure (ITT|) drawn by taking P3 = P4 = 
shows that the amplitude of the orbit increases 
with time. 

6. Stability Criteria 

Since, It is already known that to find the 
orbital parameters of the infinitesimal mass, we 
require position {x, y) as well as velocity co- 
ordinates (i, y) at different instant i.e. we have a 
four dimensional phase space (x, y, i, y). Again, if 
we take x, y and x as our three variables then the 
other one y, is obtained with the help of Jacobi 
integral equation ([7|) provided that the value of 
Jacobi constant at time i = is known. Now, 
we define plane y — Q va the remaining three di- 
mensional space (x,2/,i), then the plot x Vs x 
at every time when the trajectory intersects the 
plane y = in a part icular direction y > gives 
the PSS (|Murrav and Dcrmott 2000[). The PSS is 
a very helpful to determine the region of stability 
of the trajectory. The smooth well defined islands 
represent the trajectory is likely to be regular 
whereas the islands liberate an exact resonance in 
between the mean motions of infinitesimal mass 
and perturber. In other words, appearance of such 
islands is a characteristic of resonant motion. In 
that cases a mean motion resonance of the form 
0.1 + 0,2 : ai, where oi and 02 are integers, provides 
02 islands. Any fuzzy distribution of the points in 
the PSS indicates that the trajectory is chaotic. 
According to KAM theory, when a curve shrink 
down to a point then this point corresponds to a 
periodic orbit. In general, the regular regions in 
the PSS represents the stability regions but inside 
the regular region there is a negligible chaotic re- 
gion whereas out side of that region motion will 
always be unstable. 
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Fig. 12. — Poincare surfaces of section at different 
values of Jacobi constant: (a)C = 2.70, {b)C = 
2.80, (c)C = 2.90, (d)C = 3.00, (e)C = 
3.10, (/)C = 3.20, {g)C = 3.30 and {h)C = 3.35 



Figure (IT2|) shows the PSS at various values of 
Jacobi constant C. From frames (a) to (d), it is 
clear that the regular islands correspond to the 
stability region get increases from C — 2.70 to 
C = 3.00 while in frames (e) to (h), it seems to 
be disappear gradually from C — 3.10 to C = 
3.35. Also, the fuzzy distribution of the points 
in the PSS decrease with an increase in the value 
of C. In other words, extent of chaos decreases 
with C. The existence of chaotic regions does not 
mean that the trajectory is unbounded but they 
cover a large area of the plane (x, i). For example, 
frames (g & h) shows that the orbits can remain 
confined for a long time at a particular value of 
C. The PSS in figure (JSl), plotted for the initial 
conditions x and x lie in the interval (—0.4, 0.4) in 
which frames (a) and (b) correspond to two sets of 
parameters. Clearly, It can be seen that there are 
two islands in the plane (x, x) corresponds to each 
of the initial conditions within the range which 
indicates that orbits have mean motion resonance 
of second order. 

7. Conclusion 

We have studied different aspects of the peri- 
odic motions of infinitesimal mass in the neighbor- 
hood of equilibrium points for the mass parameter 
< /i < 5 and found the periodic orbits in addi- 
tion to some other kind of orbit like hyperbolic, 
asymptotic etc. We have seen that some of peri- 
odic orbits are bounded and have elliptic shapes 
while some are unbounded and have spiral shapes. 
We have also found that the motion of infinites- 
imal mass is affected by radiation factor as well 
as oblateness coefficients of the primaries respec- 
tively. The stability of the motion is discussed 



Fig. 13. — Poincare surfaces of section at C = 
3.10: (a) when qi = 0.75 and A2 = 0.0 (b) when 
qi = 1.0 and A2 = 0.0025 
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with the help of PSS method and found that the 
stabihty regions first increases with Jacobi con- 
stant from C = 2.70 to C = 3.00 and then it 
decreases gradually from C = 3.10 to C = 3.35. 
Thus, we conclude that the perturbation factors 
and the presence of disk play a very significant role 
for the periodic motion of the bodies in space. 
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